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O ' Abstract 

i> : 

^^ , We identify a set of quantum graphs with unique and precisely 

defined spectral properties called regular quantum graphs. Although 
chaotic in their classical limit with positive topological entropy, regu- 
lar quantum graphs are explicitly solvable. The proof is constructive: 
we present exact periodic orbit expansions for individual energy lev- 
P5 ■ els, thus obtaining an analytical solution for the spectrum of regular 

^ . quantum graphs that is complete, explicit and exact. 

^' PACS: 05.45.Mt,03.65.Sq 

A quantum graph ||l|, H, |, §, H, ||, |^, H is a network of vertices and bonds 
/\ • with a quantum particle moving along its bonds. An example of a graph with 

j^ ■ five vertices and six bonds is shown in Fig. 1. The wave function and the 

energy levels of a quantum particle on a graph are defined by the correspond- 
ing one-dimensional Schrodinger equation. Despite the apparent simplicity 
of the system, quantum graphs have proven to be a rich source of physical 
insight. From the mathematical point of view, the spectral properties of 
Schrodinger operators on graphs are highly nontrivial and have been widely 
investigated in the mathematical literature |]l], 0, 0, Q. Quantum graphs 
and networks have also been used to model various phenomena in different 
branches of physics and chemistry for more than 30 years. The most recent 
physical development appeared in a series of publications [H, P, 0, |H , where 
quantum graphs were studied in the context of quantum chaos. 



It is easy to see that the behavior of a particle on a quantum graph is 
very complex. Each time the particle encounters a vertex Vi of the graph, it 
can scatter with different probabilities in the forward or backward directions 
along any of the bonds emanating from the vertex. As a simple physical 
analogue of this system one may imagine a beam of light traveling along a 
network of optical fibers. At every joint of the fibers the light waves scatter 
in such a way that the total energy flux is conserved. 

As a result of the multiple scattering possibilities at the vertices, the 
dynamics of a classical particle on a graph is very complex and the number 
of possible periodic orbits traced by the particle increases exponentially with 
their lengths. Consequently the topological entropy of the particle is positive 
and since the phase space of the system is bounded, the dynamics of the 
particle is ergodic. The classical chaoticity notwithstanding it was shown 
that several important spectral characteristics of quantum graphs, such as 
the density of states and the spectral staircase, can be obtained exactly in 
terms of periodic orbit expansion series [H H 0- 

The "wiring" of a quantum graph, i.e. the arrangement of bonds and 
vertices, is called the topology of the quantum graph. For any given graph 
topology there exists a wide variety of possible quantum graphs. The vertices, 
for instance, may be realized as simple hubs that redistribute the quantum 
flux into various channels; or we may place 5 function potentials of vari- 
ous strengths at the vertices; or there can even be potential functions on 
the bonds. Because of this flexibihty we anticipate that the topology of a 
quantum graph alone does not uniquely specify the graph's spectral proper- 
ties. We conjecture that there are several different types of quantum graphs, 
all chaotic in their classical limit, but each exhibiting unique and precisely 
defined spectral characteristics. 

We start a rigorous classification of quantum graphs by defining regular 
quantum graphs. Although generically their classical limit is chaotic, we 
show that their spectrum is explicitly solvable analytically, state by state, 
via explicit periodic orbit expansions. This result is backed up by rigorous 
mathematical proofs [Q whose basic elements are presented below. To the 
authors' knowledge this is the first time that the spectrum of a quantum 
chaotic system is obtained exactly and explicitly. In addition regular quantum 
graphs show a spectral gap at small energy spacings and a cut-off at large 
spacings. Both the size of the gap and the location of the cut-off are computed 
analytically. 



The regular quantum graphs are different from a class of quantum graphs 
described in the hterature |^, ||, whose level spacing distributions, to a 
good approximation, exhibit features of the Gaussian orthogonal (unitary) 



ensemble |T^, |TT|. Thus there are at least two different types of quantum 
graphs with distinct spectral properties. 

The periodic orbits of graphs are defined as the periodic connected se- 
quences of bonds Bij. Denoting by k the wave number (momentum) of the 
particle, each periodic orbit contributes 

Sij = / k{x) dx, (1) 

J Bij 

to the total action of a path traced by the particle. It turns out that the 
information contained in the totality of all the possible classical periodic 
orbits often allows one to reconstruct exactly certain quantities of a purely 
quantum nature. For example, according to recent results 0, |^, |^, |^, the 
exact periodic orbit expansion for the density of states can be written as 
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p(fc)^5:<5(fc-fc,)=p(fc) + -Re^r,^A;:e-^^('=). (2) 

j=l ^ p u=l 

Here p{k) is the average density of states, u is the repetition index, Tp = 
dSp{k)/dk, and Sp, Ap are correspondingly the action and the weight factor 
of the prime periodic orbit labeled by p. We assume in what follows that the 



system is scaling g, |1^, |13|, |Tj, ^, |T6|, |T3- This means that the actions of 



the periodic orbits are proportional to the wave number, 

Spik) = S'pk, (3) 

where 5°, a constant, is the reduced action. In this case Tp = Sp and Ap 
are /c-independent constants. We define the total reduced action So of the 
graph, 5*0 = {J kdx)/k, where the integral is extended over all of the bonds 
of the graph. The scaling assumption is not an artificial restriction. It 
occurs, for instance, in atomic physics where scaled spectroscopy is now a 
common experimental technique |T^. In addition, scaling quantum systems 



of this kind are the analogues of certain electromagnetic ray-splitting systems, 
fiat metal cavities partially filled with a dielectric substance such as Teflon 



The formal description of a quantum graph system proceeds as follows 
[^, H, 1^. On a bond B^ connecting the vertices Vi and Vj, the wave function 
of a quantum particle is defined by the one- dimensional Schrodinger equa- 
tion which may include potentials on the bonds |T^. At every vertex Vi the 



wave function satisfies the usual boundary conditions of continuity and flux 
conservation. The consistency of the boundary conditions at every vertex 
of the graph naturally yields the global quantization conditions that deter- 
mine the momentum eigenvalues kn- A simple and elegant method based on 
the scattering quantization approach was presented in [^, H, 0], where the 
quantization condition is given in the form 

det[l - S{k)] = 0. (4) 

Here S{k) is the scattering matrix [^, written explicitly in terms of the 
connectivity matrix [Q of the graph. It can be shown that the modulus of 
the complex function (^ is a trigonometric function of the form 

cos(5oA; + 7r7) - <!>{k) = 0, (5) 

where 

$(A;) =J2ai cos{Sik + 7r7i). (6) 

i 

In the scaling case the Oj are constants and Si < Sq are certain combinations 
of the classical actions (^. In general the functions 7 and ■ji are bounded 
and tend to constant values for k —>■ 00. In the scaling case they are k- 
independent constants. 

In order to proceed we define regular quantum graphs. A regular quantum 
graph fulfils the condition 



a 



E I «'l < 1- (7) 



It is convenient for the following discussion to assume a > 0, i.e. we exclude 
trivial graphs with a = 0. They are regular quantum graphs whose spectrum 
can be obtained trivially. For regular quantum graphs (j^) can be solved 
formally for the momentum eigenvalues kn- We obtain the following implicit 
solution for the roots of this quasi periodic function 

, ^ r , 1,1 rarccos($„), for n + /i even ,„. 

5^0 So I TT — arccos (<!>„), tor n + /x odd, ^ ' 



where $„ = ^{kn) and fi & Z, a fixed integer, is to be chosen such that ki 
is the first nonnegative solution of @. Because of (|^), the boundedness of 
the trigonometric functions in (^) and the properties of the arccos function, 
the second term of (H) is bounded away from and vr/S'o and assumes only 
values between u and n/So — u, where < u = arccos(a)/S'o < n /2Sq. Thus, 
for regular graphs the points 

A:n = ^(n + /x-7 + l), n = l,2,... (9) 

'-'0 

are not roots of (p. They serve as separators between root number n and 
root number n + 1. In fact (^ implies even more: the existence of finite- 
width "root-free zones" Rn = (kn — u,kn + u) surrounding every separating 
point fcn, where no roots of (|) can be found. Thus, roots of (H) can only 
be found in the "root zones" Z„ = [fc„_i + «,/;;„ — «], subsets of the root 
intervals /„ = [/;;«- 1, A^n]- Since Sq is the largest action in (||) and (^, it can 
be shown that kn is the only root in Zn- Therefore, in summary, there is 
exactly one root kn inside of Zn C /„, and this root is bounded away from 
the separating points kn-i and kn by a finite amount u. 

The spectral properties of regular graphs discussed above allow us to 
draw several important conclusions, (i) Since there is exactly one root kn 
of @) in Ini this proves rigorously that the number of roots of @) smaller 
than k grows like N{k) ~ Sok/ir (Weyl's law), while the second term of 
(H) represents the fluctuations around the average, (ii) The existence of the 
root-free zones Rn gives rise to a spectral gap of finite width g = 2u in 
the nearest neighbor spacing distribution of regular quantum graphs, i.e. no 
level spacings smaller than g are ever found. This, for instance, precludes 
the possibility of degenerate eigenvalues, (iii) The existence of the separating 
points (^ together with the root-free zones imply the existence of a spectral 
cut at c = 2(7r/S'o — u), i.e. no level spacings greater than c are ever found. 
(iv) The existence of the separating points @ and the root-free zones Rn are 
the key for obtaining an explicit and exact periodic orbit expansion for every 
root of @. Multiplying both sides of (0) by k and integrating from kn-i to 
kn we obtain 
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kn = kn- -^ -Re 2^2^Ap 

^ p u=l 
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where Up = TrSp/So- Since all the quantities on the right-hand side of (p!0|) 
are known, this formula provides an explicit representation of the roots kn 
of the spectral equation (^ in terms of the geometric characteristics of the 
graph. To our knowledge, this is the first time that the energy levels of a 
chaotic system are expressed explicitly in terms of a periodic orbit expansion. 
Previously, explicit formulae for individual energy levels were known only for 
integrable systems. In the context of periodic orbit theory, the energy levels of 
integrable systems are given by the Einstein-Brillouin-Keller (EBK) formula 
[ |T^ . However, apart from a few exceptional cases |T^ EBK quantization is 
only of semiclassical accuracy. 

For a generic chaotic system the energy levels can only be obtained in- 
directly as the singularities of the periodic orbit expansion of the density of 
states (Gutzwiller's formula ||18[), an implicit method which, in general, is 
only of semiclassical accuracy. Formula ([TI1|), on the other hand, shows that 
for regular quantum graphs every quantum level can be obtained individually, 
explicitly and exactly in terms of classical parameters. 

In order to demonstrate that the class of regular quantum graphs is not 
empty we present an explicit example: the one- dimensional step potential 



shown in Fig. 2 (a) |20|. With 

V^o = AE (11) 

we turn it into a scaling system. Physically this potential is realized for 
a rectangular microwave cavity partially loaded with a dielectric substance 



1^, |T^, |T5| . The scaling step potential is equivalent to the scaling three- vertex 
chain graph shown in Fig. 2 (b), with two bonds Li = b and L2 = /3(1 — b), 
(3 = Vl — A. In this case the spectral equation is given by 

det [1 - S{k)] = sm{Lk) - r sin [(Li - L2) k] = 0, (12) 

where L = Li + L2, and r = {1 — (3) / {1 + f3) is the reflection coefficient at 
the vertex V2 between the two bonds. The solution of (0) is 

TTTl 1 

fc„ = — + (-1)"-^ arcsin [r sin((Li - L2) A;„)] , (13) 

where the second term, which gives the fluctuations around the average, is 
bounded by 

1 IT 

l^^nl < y arcsin(r) < — . (14) 
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Therefore, for the three- vertex chain graph with r < 1, the separating points 
(H) are given by 

kn = T(^+l)- (15) 



They mark the borders of the root intervals /„ for all possible values of the 
bond lengths. For a given prime periodic orbit p, such as the ones shown in 
Fig. 2, the coefficients Ap are given by P, P^ 



Aj, = (-l)x(p)^-(p)(i _ rY^p)/2^ (16) 

where r is the reflection coefficient at the middle vertex and a{p) and r(p) are 
correspondingly the number of reflections off it, and transmissions through 
it. Since the reflection coefficient may be positive or negative depending 
on whether the particle scatters from the right or from the left, the factor 
(_X)x(p) ig needed to keep track of how many times it appears with a minus 
sign. Moreover x keeps track of how many times the particle scatters off 
the walls of the potential, since each scattering event from the walls gives 
rise to a minus sign in the wave function. Using the explicit form of the 
spectral determinant (0), one obtains the explicit form of the eigenvalues 
for a quantum particle in the step potential in terms of classical (Newtonian 
and non- Newtonian |]13[) periodic orbits. 

In order to illustrate the validity of (p!0D we computed ki, kio and kioo 
including periodic orbits that experience up to five scattering events. We 
chose b = 0.3 (see Fig. 2) and A = 1/2 (see ([11])). We obtain fc^ = 
4.1161, fciJ = 39.2866 and k[% = 394.9477. This can be compared with 
the exact k values given by ki = 4.107149, kio = 39.305209 and /cioo = 
394.964713. In order to illustrate convergence, we also computed ki, kio and 
fcioo including all periodic orbits that experience up to 20 scattering events. 
This amounts to including more than 100,000 periodic orbits and results 
in fcf°^ = 4.105130, k^^^ = 39.305212 and fcioo = 394.964555. This also 
indicates that the convergence of (p!0| ) is not destroyed by keeping more and 
longer periodic orbits. 

Additional examples of regular graphs are provided by the scaling "Man- 
hattan potentials" . These potentials are obvious generalizations of the step 
potential shown in Fig. 2 (a) to piecewise constant potentials where the po- 
tential heights scale with the energy. Furthermore we checked explicitly that 
linear chain graphs with scaling 6 function potentials at the vertices provide 
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more examples of regular quantum graphs. In this case the strengths of the 
6 function potentials scale linearly with the momentum. 

The key physical feature of these one- dimensional quantum systems, which 
permits the exact periodic orbit expansion ([T0[) for the eigenvalues, is the 
rigidity of their spectra. For integrable systems, spectral rigidity is due to 
the "geometrical rigidity" of the periodic orbits, confined to integrable tori. 
In the case of quantum graphs, the geometrical structure of the periodic orbits 
is much more complicated. The complexity of the expansion (p!OD compared 
to the EBK formula reflects the geometrical complexity of the periodic orbits. 

This paper profitted from a fruitful exchange of E-mails with Uzy Smilan- 
sky. Y.D. and R.B. gratefully acknowledge financial support by NSF grants 
No. PHY-9900730 and PHY-9984075; Y.D. and R.J by NSF grant No. PHY- 
9900746. 
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Figure Captions 

Fig. 1: A generic (quantum) graph with five vertices and six bonds. 

Fig. 2: (a) Simple step potential, a basic problem in one-dimensional quan- 
tum mechanics. Also shown are examples of Newtonian ("N") and non- 
Newtonian ( "NN" ) periodic orbits used in the periodic orbit expansion of its 
energy eigenvalues (see Text), (b) Three- vertex chain graph corresponding 
to the step potential above. 
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